Robustness of dynamical decoupling sequences 
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Active protection of quantum states is an essential prerequisite for the implementation of quantum 
computing. Dynamical decoupling (DD) is a promising approach that applies sequences of control 
pulses to the system in order to reduce the adverse effect of system-environment interactions. Since 
every hardware device has finite precision, the errors of the DD control pulses can themselves 
destroy the stored information rather than protect it. We experimentally compare the performance 
of different DD sequences in the presence of an environment that was chosen such that all relevant 
DD sequences can equally suppress its effect on the system. Under these conditions, the remaining 
decay of the qubits under DD allows us to compare very precisely the robustness of the different 
DD sequences with respect to imperfections of the control pulses. 



I. INTRODUCTION 

Quantum computers can execute certain tasks more 
efficiently than classical computers by processing infor- 
mation according to the laws of quantum mechanics. In 
analogy to a classical bit, which can assume the values 
or 1, quantum mechanical two-level systems like a spin- 
1/2 can be used as quantum bits by identifying their 
eigenstates with these values, e.g. |0) for spin up and 
|1) for spin down. In quantum information processing 
and quantum memory applications, it is very important 
to keep the information isolated from the environment: 
uncontrolled interactions with the environment tend to 
degrade the quantum information. This environment- 
induced loss of quantum information is called 'decoher- 
ence' [1-3]. 

If one is able to control the system in such a way to 
reduce the detrimental effect of the system-environment 
(SE) interaction, one can preserve the quantum state for 
a longer time. This way of fighting decoherence by ap- 
plying fast and strong pulses has been termed dynami- 
cal decoupling (DD)[ - ]. The main attraction of DD is 
that it requires few additional resources. In contrast to 
quantum error correction [10-12], e.g., it does not require 
additional qubits. DD can be traced back to Hahn's 'spin 
echo' experiment, where a refocusing pulse induces a time 
reversal of the SE interaction of nuclear spins [13]. This 
increases the decay time or decoherence time of the stored 
information in the qubit. The technique has evolved sig- 
nificantly since then and its efficiency was studied and 
demonstrated in many different systems [14-30]. 

It has been shown that the type as well as the spectral 
density of the SE interaction play a significant role for 
finding the optimal DD sequences [7, 25, 31-38]. More- 
over, unavoidable errors in the control pulses are also an 
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important source of decoherence [18, 24, 39-49]. Thus 
optimal DD sequences must be able to reduce the ef- 
fective SE interaction while compensating the effects of 
non-ideal control fields [22, 24, 43, 46-51]. Reference [48] 
is a recent review of this subject. 

In this article we compare the performance of differ- 
ent DD sequences in a system where pulse errors are 
the dominant source of decoherence. To the best of our 
knowledge, this situation has not yet been investigated: 
all previous works were done in systems where the pulse 
errors only become significant when most of the effects 
of the SE interaction have been eliminated. To this end, 
we prepare a system where the spectral density of the 
SE coupling has two main contributions. One source of 
noise is almost static and can therefore be refocused by all 
tested DD sequences. The other contribution is a rapidly 
fluctuating noise, whose correlation time is much shorter 
than the time required for an inversion of the spins. This 
type of noise cannot be refocused by any DD sequence. 
Therefore the main difference between the performance 
of the DD sequences is their susceptibility to pulse errors. 

In a previous work [18, 24], we compared the perfor- 
mance of different DD sequences and in particular we 
found that the 'Carr-Purcell-Meiboom-GiU' (CPMG) 
sequence [52, 53] performed particularly well for specific 
initial states. In this case the decoherence time was one 
order of magnitude longer than for robust sequences that 
reduced decoherence symmetrically with respect to arbi- 
trary initial states. We found that this difference arose 
because the system qubit, a 13 C spin, interacted with 
neighboring 13 C spins. While the CPMG sequence was 
able to reduce the effect of 13 C- 13 C couplings, the robust 
sequences were not [36]. Since the tested DD sequences 
were not designed for eliminating the effect of homonu- 
clear couplings, the longer decoherence times for CPMG 
applied to certain initial conditions cannot be taken as 
a measure of its performance. In this work we therefore 
use a system that does not exhibit homonuclear interac- 
tions. This allows us to show that the robust sequences 
can also achieve the optimal decoherence time observed 
under CPMG, and this performance is independent of the 
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initial condition. 

The article is organized as follows. In Sec. II we define 
the system, in Sec. Ill we introduce the DD sequences to 
be compared and in Sec. IV we compare their robustness 
against pulse errors effects. In the last section we draw 
some conclusions. 



II. THE SYSTEM 

The experimental system is an ensemble of non- 
interacting spins 1/2. They consist of the protons of a 
water sample to which we added 5 mg/100 ml CUSO4 to 
reduce the T\ relaxation time to 287 ms. This results in 
faster repetition times and shorter overall duration of the 
experiments. The sample was placed in a static magnetic 
field along the z-direction and its Hamiltonian is 
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where uj s is the Zeeman frequency and S z is the sys- 
tem spin operator along the z axis. The inhomogeneities 
of the static field correspond to a static perturbation, 
and molecular motion makes this perturbation time- 
dependent on a time scale that is slow compared with the 
delays between the DD pulses used in our experiments. 
This makes it possible to refocus this perturbation very 
effectively. 

The second major source of noise is the fluctuating 
dipole-dipole interaction, whose correlation time is the 
molecular reorientation time (« 35 fs) - much faster than 
any conceivable control fields for nuclear spins and there- 
fore not amenable to DD. On the other hand, these fluc- 
tuations are so fast that their average effect on the system 
is relatively small [54]. 

Experiments were performed on a home-built NMR 
spectrometer with a 'H resonance frequency of 360 MHz. 
The radio frequency field strength was 2ir ■ 13.3 kHz, 
which corresponds to a n~ pulse duration of 37.5 us. An 
initial state oc I x or I y was prepared by rotating the I z 
equilibrium state with a resonant 7r/2 pulse. The free 
evolution decay of the transversal magnetization of our 
system has a decay time of 2.9 ms (free induction decay). 
A simple Hahn-echo sequence [13] increases this time to 
106 ms, as shown in Fig. 1. 



III. DYNAMICAL DECOUPLING SEQUENCES 



Hahn echo decay 
■Exponential fit (HE) 
Free evolution 
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Figure 1. (Color online) Decay of the magnetization of the 
1 H spin system under free evolution and in a Hahn spin-echo 
(HE) sequence. 
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DD sequences consist of repetitive trains of 7r-pulses. 
The delays between the pulses and their phases are im- 
portant parameters for improving the performance of the 
DD sequences [4-6, 50, 52, 53, 55, 56]. In particular the 
relative phases, which correspond to the directions of the 
rotation axes, are important for making the sequences 
robust against pulse imperfections and unwanted envi- 
ronmental interactions [18, 22, 24]. 

Fig. 2 gives an overview of the sequences that we ex- 
amined for this work. It shows a single cycle for each 



Figure 2. (Color online) Dynamical decoupling pulse se- 
quences tested in this work. 



sequence, which is repeated as often as required, r is 
the delay between the pulses. The Carr - Purcell (CP) 
sequence [52] and the version by Meiboom and Gill [53], 
known as CPMG, use the same sequence of refocusing 
pulses; they only differ with respect to the state to which 
they are applied. In the case of the CP sequence, the 
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initial state is perpendicular to the rotation axis of the 
inversion pulses, in the CPMG version, it is parallel. Er- 
rors in the flip angles destroy the perpendicular com- 
ponent, but they leave the longitudinal component un- 
scathed [18, 53]. The sequence XY4 was introduced by 
Maudsley [55] and it reduces the effect of pulse imperfec- 
tions for arbitrary initial states [18, 55, 56]. It consists 
of four pulses with phases x — y — x — y (Fig. 2). An 
asymmetric version of the XY-4 sequence was introduced 
by Viola et al. [ ], which we designate XY4(a). The 
XY8-sequences are symmetrized versions of the XY4 se- 
quences [4G, 56]. Two DD sequences that are particularly 
robust against flip-angle and resonance offset errors are 
the KDDa; and KDD^ sequences [24, 47, 48] . They were 
designed by combining the rotation pattern of the XY4 
sequence with that of a robust composite pulse [57]. 

In earlier works using these sequences, the conditions 
were chosen such that the dominant perturbation was the 
environmental noise [18, 24, 44, 46, 47]. In this work, we 
focus on a system that allow us to make a comparison 
between these sequences in a regime where all sequences 
perform equally well at eliminating the environmental 
noise and any differences in their performance can be 
attributed directly to their robustness, i.e. to their effi- 
ciency in suppressing the effect of pulse imperfections. 



IV. ROBUSTNESS COMPARISON 

To compare the sensitivity of the sequences to pulse 
imperfections, we prepared two orthogonal initial states 
I x and I y and then measured their decay as a function of 
time under the application of the different DD sequences 
described in the previous section. Figure 3 shows the echo 
train of a CPMG sequence. From this data, we extracted 
the signal (I x magnetization in this example) at the end 
of each DD cycle (marked by blue squares in the figure). 
The decay of the echoes was mostly exponential, with 
some exceptions discussed below. 

The extracted echoes (blue squares in Fig. 3) were fit- 
ted with an exponential function to obtain the decay time 
of the magnetization. Experiments were repeated and the 
decay times plotted as a function of the delay between 
pulses. 

Figure 4 compares the decay times of different DD se- 
quences as a function of the delay between pulses. For 
the CPMG sequence, we present the decay of the I x and 
I y magnetization separately, marked as CP and CPMG, 
respectively. For the other sequences, whose performance 
is quite symmetric with respect to the initial condition, 
we present the decay times averaged over the two initial 
conditions. For long delays between the pulses, the ob- 
served decay times reach a limiting value of ~ 276 ms, 
irrespective of the sequence and the initial condition, and 
very close to the measured value of T\ . This is a verifica- 
tion of the assumption that all sequences can effectively 
decouple the slowly fluctuating environment. 

For shorter pulse delays (i.e. more pulses in a given 
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Figure 3. (Color online) Time evolution of the spin-system 
magnetization under the application of a CPMG sequence. 
The black solid line shows the evolution of the magnetization 
and the blue squares mark the echo amplitude at the end of 
a CPMG cycle. We use the echo maxima for measuring the 
CPMG decay time. The cycle time was r c = 32 ms. 
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Figure 4. (Color online) Average decay times as a function of 
the delay r between pulses for different DD sequences . 

time interval), the signal decays more rapidly. This is 
most prominent for the CP sequence. In this situation, 
pulse imperfections add coherently and generate a rapid 
loss of magnetization [18] . 

As the pulse delays become shorter than 0.5 ms, which 
corresponds to 864 pulses during the 0.5 s measurement 
time, the other sequences also start to generate shorter 
decay times, and their decays become nonexponential. 
Figure 5 shows a representative example of such a signal. 
It can be fitted with a double exponential, 

s{t) = ae- t ' T i +6e-*/ r »' (2) 

with two decay times and T|. 

Figure 6 shows the decay times fitted with eq. (2) for 
different average delays between pulses. For the fast de- 
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Figure 5. (Color online) Normalized spin magnetization as a 
function of the evolution time for short delays (t = 100 fis) 
between the pulses for an XY8(s) sequence. Pulse errors dom- 
inate here, inducing a multi-exponential decay. The red solid 
line is a fit to eq. (2). 
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Figure 6. (Color online) Fitted decay times for the double- 
exponential decay 



resulting per pulse of the sequence. For this evaluation, 
we consider only the short-time component described by 

Since the pulse error is the dominant source of decay, 
we quantify its effect by measuring the fractional decay 
of the magnetization per pulse. The pulses are the same 
for all the DD sequences, but their effect, averaged over 
full cycles, shows how well the sequence is able to cancel 
the imperfections of the individual pulses. 



CD 
U3 
=J 
Q_ 
i— 
CD 
Q_ 



CD 

CD 
O) 
CO 

V- 

CD 
> 
< 



10" 



10" 



10" 




XY4(a) 
XY4(s) 
XY8(a) 
XY8(s) 
KDD 

X 

KDD 



CPMG 



10 20 30 40 50 60 70 80 

Number of pulses 

Figure 7. (Color online) Average error per pulse for different 
DD sequences with delay r —lOOfis. 

Figure 7 shows the average decay per pulse for the dif- 
ferent sequences, plotted against the number of pulses. 
For these data, the interpulse delay was t = 100 (is. The 
most conspicuous feature is that CP performs very badly 
and CPMG very well. The compensated sequences lie be- 
tween these two extremes, and we find that the higher or- 
der sequences (XY8, KDD perform better than the lower 
order sequences (XY4). For unknown initial conditions, 
KDD shows the best performance. Under the present 
conditions, sequences that differ only with respect to time 
reversal symmetry perform quite similarly, in contrast to 
other cases discussed earlier [46]. 



cay times (T%), which are represented by empty sym- 
bols in Fig. 6, the performance of all DD sequences is 
quite similar. For the slow component (T|, represented 
by filled symbols), XY8, CPMG, and KDD I perform bet- 
ter than XY4. 

As we stated before, the decay time is dominated by 
the effect of pulse errors for short delays. The resulting 
average Hamiltonian projects the magnetization onto its 
eigenbase - this results in the fast decay component. Af- 
ter this projection, the remaining magnetization, which 
is not significantly affected by the pulse imperfections, 
decays on a slower timescale, which is dominated by the 
environment. 

For a quantitative comparison of the different pulse 
sequences, we calculate the average magnetization decay 



V. CONCLUSION 

We have have tested the robustness of different DD 
sequences by comparing them in an environment that in- 
teracts with the spins in such a way that the decoherence 
time under the application of DD sequences with ideal 
pulses is independent of the delay between the pulses. 
This allowed us to study the robustness of the different 
DD sequences by isolating the effects of the pulse errors. 
We found that the decoherence time of the most robust 
sequences, the KDD family, is the longest for arbitrary 
initial states. This is consistent with the measured error 
per pulse averaged over many cycles of a DD sequence, 
where the KDD sequences have the lowest effective er- 
ror. In the regime studied, the time symmetrization on 
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the cycles does not play a significant role for reducing for helpful discussion. This work was supported by the 

decoherence, since only the phases of the pulses are im- DFG through grant Su 192/24-1. 
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